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If the speed of sound were vastly larger in the early Universe a near scale-invariant spectrum of 
density fluctuations could have been produced even if the Universe did not submit to conventional 
solutions to the horizon problem. We examine how the mechanism works, presenting full mathemat- 
ical solutions and their heuristics. We then discuss several concrete models based on scalar fields 
and hydrodynamical matter which realize this mechanism, but stress that the proposed mechanism 
is more fundamental and general. 
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PACS numbers: 0000000 

1 . Introduction The fact that the large scales we ob- 
serve today were "out of touch" in the early Universe 
is one of the greatest annoyances of standard cosmology. 
This "horizon" problem prevents a causal explanation for 
the observed features of our Universe, which have to be 
added on "by hand" as initial conditions. Nowhere is this 
more unsatisfactory than in relation to the primordial 
density fluctuations that seed the structures we observe 
today. We have measured these structures with tremen- 
dous accuracy; yet the primordial fluctuations cannot be 
explained in the standard scenario and remain God-given 
initial conditions. This weakness has motivated many re- 
visions to the standard Big Bang picture, notably infla- 
tion [j. 

Any explanation for the initial spectrum of fluctuations 
has to begin with a mechanism for causally connecting 
vast scales at the cradle. But this is barely the beginning: 
one then has to suggest a physical mechanism that would 
render these scales homogeneous but "not quite" . In the 
detail lies the hurdle: the inhomogeneities must be near 
scale-invariant, i.e. look approximately the same on all 
scales, and have a well-defined amplitude, of the order 
of a part in 100,000. Herein lies the challenge for any 
structure formation scenario. 

Solutions to the horizon problem rely on either accel- 
erated expansion a contracting phase followed by a 
bounce Q, or a loitering early phase An alterna- 
tive was supplied by varying speed of light (VSL) theo- 
ries 0, : the idea that the early Universe operated 
under a much larger maximal speed, causally connecting 
everything. There are several implementations of this 
idea and structure formation has been considered in a 
few Q. But in spite of these valiant efforts it is fair to 
say that VSL scenarios have so far failed to explain the 
cosmic structures. 

In this paper we use a varying speed of sound as a 
proxy for VSL. In theories containing two metrics the 
speed of massless matter particles and the speed of grav- 
ity are different. This defines a frame where gravity is 
unaffected, but "light" travels much faster. We propose 
a new mechanism for producing scale-invariant fluctua- 
tions, based on a varying speed of sound. The idea mim- 
ics that of a varying speed of light. If the sound hori- 
zon shrinks, modes that start off oscillating eventually 



freeze-out. The universe is initially homogeneous apart 
from small quantum or thermal fluctuations. These are 
left imprinted "outside the horizon" after the speed of 
sound c s has decreased suitably. 

A varying speed of sound has been considered in 
other guises before. We nod at K-essence models, 
based on (scalar) field theories with non-standard kinetic 
terms [20]. These have been used as inflationary and 
quintessence models, but can be adapted to implement 
the varying speed of sound mechanism we advocate. The 
cuscuton model of [13l |14| provides a related framework. 
We can also bypass scalar fields, and simply regard the 
speed of sound c s and the equation of state w as free phe- 
nomenological parameters [l6|, [l7| • Even though VSL is 
our leading motivation here, we stress that many other 
approaches may be linked with the conclusions in this pa- 
per. A VSL implementation, however, may be required if 
additionally one wants to solve the flatness, entropy and 
the other problems of Big Bang cosmology. 

2. The mechanism We shall first illustrate the mech- 
anism in its simplest realization. Suppose that gravity 
remains unmodified and that we restrict ourselves to ex- 
panding Universes with w > —1/3. This is to avoid con- 
fusion with inflation and ekpyrotic scenarios (although 
constructive alliances should be investigated). The new 
ingredient, now, is the assumption that the speed of 
sound is density dependent and diverges with conformal 
time like c s oc r\~ a (with a > 0; note that r\ is posi- 
tive and increases from zero). A reparameterization in 
terms of the density will be examined later (with a strik- 
ing result) but this is the most suitable expression for 
a mathematical solution. Concrete models realizing this 
set up will be presented but we do not want to wed what 
follows to any one of them. 

Whether we employ a fluid or a scalar field the density 
fluctuations are described by a modified harmonic oscil- 
lator equation. This can be written in terms of variables 
related to the Newtonian potential $ or the curvature 
perturbation^, the so-called "u" and "v" . The equation 
for v is [H, El, HU] : 
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where z cx a/c s . Thus c s appears in two places: in the k 2 
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term responsible for sub-horizon oscillations and in the 
variable mass term, in z" / z. The equivalent equation for 
u has a variable mass term that does not contain c s . 

As with inflation this equation can be exactly solved 
with Bessel functions, but we'll first examine a WKB so- 
lution in order to establish the initial conditions. With 
a > f modes start inside the sound horizon but even- 
tually leave it (the horizon scale is set by c s krj ~ f, 
i.e. k oc n 1 ^ ). We can therefore initially ignore the 
term in z" / z. Even though the frequency to = c s k is 
changing, the WKB condition uj' <C uJ 2 translates into 
c s kr\ 3> \a\, which is always satisfied early on. This 
amounts to preserving the adiabatic invariants, i.e. the 
number of quanta in a given mode is kept fixed but not 
the total energy in that mode, which changes propor- 
tionally to the frequency. The appropriately normalized 
WKB solution is therefore 



Jk J c s d-q e -pc e kq 



(2) 



where (3 = l/(a — 1) > and in the last ~ we neglected a 
phase. We can also consider scenarios where c s is initially 
constant, then drops like a power law. The standard 
boundary condition may then be imposed in the constant 
c s phase and propagated using the WKB solution, with 
the same result . 

As with the equivalent calculation in inflation, Eqn. ([!]) 
can be transformed into a Bessel equation, with solutions: 



v = ^/f3n(AJ u (f3c s kn) + BJ_ u (0c s kr))) 



(3) 



The order v is given by v = (3 (^a 
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— t- — - — i and A and B 

are k- independent numbers of order 1, so that the bound- 
ary condition {2J is satisfied. The spectrum left outside 
the horizon is now easy to find. Since c s r\ is a decreas- 
ing function of time, the negative order solution is the 
growing mode, so that asymptotically we have: 



(4) 
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Since the curvature fluctuation is related to v by £ = v/z, 
its scale- invariance (fc 3 £ 2 = const) requires v = 3/2, i.e. 



a = ao = 6 
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If we rephrase this requirement by writing c s in terms 
of the density p we conclude, interestingly, that c s oc p 
for all w. The spectrum can also be made red or blue 
depending on whether a < qq or a > qq, specifically 



n s - 1 = 0(a - oq) 



(6) 



We note that an infinitely fast transition (a>l) implies 
ns = 2 for all w. All of these considerations depend on 
the sub-horizon normalization, here chosen to match a 
vacuum quantum state. 



We can also work out the fluctuations' amplitude for 
near scale-invariant spectra. Considering the curvature 
£, which "freezes-in" (i.e., is time independent outside 
the horizon, even with a variable c s ), it is found [2~H | 
after straightforward algebra: 
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Those acquainted with this expression (say, from infla- 
tion) will find here a good explanation for why c s oc p 
leads to scale-invariance, even without inflation. If we 
refine the c s law to c s = cq(1 + p/p*), (where c s <=a cq at 
low-energy and is the density that triggers its diver- 
gence) we find: 
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This forces the energy scale of the varying speed of sound 
phenomenon to be a couple of orders of magnitude below 
the Planck scale, unless w 3> 1 in which case it can be 
significantly lower. 

This is not to say that this scenario is protected from 
the super-Planckian "problem" (cf. [27j for a similar af- 
fliction). Indeed, choosing w = 1/3 to fix ideas, and 
taking into account the amplitude ©, straightforward 
algebra shows that the current Hubble scale left the 
sound horizon when H^ 1 = H~ 1 {p)c s {p)ao/a{p), trans- 
lating into p/Mpi - I0 24 , or E - I0 6 A/ P; . This can 
only be evaded with a very large w > 1, leading to 
E+ < I0~ 10 A/pi so that the current Hubble scale does 
freeze- in for E < Mpi; however for the normalization to 
be correct one would then need w > I0 30 . We would not 
discount more standard scenarios (with w of order one), 
as this "problem" might be quite fictitious: it could well 
be that there is never a Planckian quantum gravity phase. 

Perhaps more interesting is the possibility that the ini- 
tial conditions are "thermal" rather than quantum vac- 
uum fluctuations as assumed above [25| . The sub- horizon 
modes then have a spectrum obtained from the above 
by a multiplicative factor of T c /k (where T c = Ta/c s 
is a constant). This reflects the Rayleigh- Jeans limit of 
the thermal occupation number n(k) and implies a sub- 
horizon white-noise spectrum. This then propagates to 
the final result (Eqns. [3]or[S]) leading to the prediction 
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Scale-invariance is now ensured by a very fast phase tran- 
sition (a > cio, i.e. something close to a step- function). 
The spectrum is always red but can also be made arbi- 
trarily close to flat. The amplitude, if departures from 
scale-invariance are small, can be found following some 
simple algebra, keeping track of all significant multiplica- 
tive constants, and is: A 2 = a{T*/M Pl f - I0~ 10 . Thus 
these models do not suffer from a trans-Planckian prob- 
lem. With a ~ 100 we have Pi,/M pl ~ I0 -16 and the cur- 
rent horizon scale leaves the horizon at a density barely 
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an order of magnitude higher than this. If anything these 
scenarios may have the opposite problem: they push 
very low if the transition is extremely steep. For example 
with a ~ 10 50 the appropriate normalization would re- 
quire ~ lGev, with the current horizon scale leaving 
the horizon at around the same energy scale. This might 
bring these scenarios within the reach of direct experi- 
mental test. Perhaps the study of high energy thermal 
plasmas [22[ could provide a direct measurement of c s , 
placing the model under laboratory test. 

To conclude our presentation of the mechanism we note 
that in general one might have both thermal and quan- 
tum contributions. In the absence of the usual inflation- 
ary justification for a vacuum state the thermal argument 
may make more sense. Also notice that in these scenarios 
the curvature ( and the potential $ both freeze outside 
the horizon and have the same spectrum and normaliza- 
tion (within factors of order 1). This is to be contrasted 
with the pathologies known to plague some scenarios [23| . 

3. Models So far we have concentrated on investigating 
a new mechanism for producing scale-invariance and on 
its phenomenology. We specifically avoided wedding it to 
any "model" (although this has to be done to solve the 
flatness and other Big Bang problems). However we'll 
happily exhibit some possible models realizing this mech- 
anism. These should be merely seen as proof of concept. 

For example, one can appeal to a scalar field 4> endowed 
with non-trivial kinetic terms [l(| EH E3, H3| , a suitable 
potential to recreate constant w (i.e. a scaling solution), 
and a varying c s . The Lagrangian for such models has the 
form Co = K(X) - V(<j>) (with X = ^d^d^cj)). From 
the stress energy tensor we find that p = K — V and 
p = 2XK x — K + V, i.e. w depends on both K and V. 
In contrast the speed of sound is given by: 
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(10) 



i.e. it depends only on K. For any power-law function 
K(X) we can always find a power-law potential V(<f) that 
ensures a constant w (scaling) solution. An example is: 



Co = V\X\X±-m 2 cj, 2 



(11) 



for which one can find a consistency relation between 
parameters X* and m, as well as initial conditions so 
that constant w solutions, for any w, can be generated. 
As is well known this model has an infinite speed of sound 
c s (cf. Ean.lTO]): it's the so-called cuscuton field [HE1- 
But if we now add to Co a new term C\ = X(X/X ir ) n ^ 1 
with n < 1/2 this term will be sub-dominant when p S> 
X*, so that the initial, constant w solution is still valid 
in this regime. However c s will no longer be infinite, 
and using Eqn. (fT0|) we find c 2 s cx X~ n cx p x ~ 2n . Thus 
for n — —1/2 we have c s cx p, as required for scale- 
invariance. Under thermal initial conditions the more 
extreme a — > — oo limit would have to be considered. 
Note that the theory has to contain a cut-off at low X 
ensuring that the kinetic term becomes linear in X at low 
X (or else the field should decay into normal matter). 



This is only an illustration, but it highlights the limi- 
tations of standard thermodynamical arguments relating 
c s and w. For hydrodynamical matter w — pj p whereas 
c 2 = Sp/Sp. If the 6 in this expression is an adiabatic 
partial derivative and the background also evolves adia- 
batically [H,EH we must have w = — 3^(1 + w)(c 2 — w), 
which results from w tP — {c 2 — w) / p, itself a simple rear- 
rangement of: 



dp 



(12) 



These relations contradict the assumptions in our per- 
turbation calculation (w cannot be constant), however 
they are blatantly violated by scalar fields endowed with 
a potential, as in the example just presented, or more 
prosaically in the cases of inflation or quintessence (for 
which c s = 1 but w can be anything). 

Once this is noted there is no reason not to consider 
general forms of matter with c s and w regarded as in- 
dependent variables [TtJ ■ Even with thermalized matter 
the evolution may be non-adiabatic [l6j either for the 
fluctuations, the background, or both, thereby violating 
(fT2")l . If the background evolves adiabatically there must 
be entropy production pT = 5p — c 2 d Sp in order to have a 
constant w and a varying c s . Specifically, the conditions 
of our calculation are met if T = (p/p+^S, which can be 
seen as the appropriate "non-adiabatic condition" . 

It is also possible that a more complicated solution to 
(JXJ) , allowing for a varying w satisfying (fT2")) and using 
the proposed varying c s mechanism, still finds a niche 
for scale-invariance. This possibility is currently un- 
der consideration, for example in the context of adap- 
tations of the Chaplygin gas (2(| (or generalizations of 
the work in [2^, Ht|). The mathematics of such mod- 
els is considerably less straightforward, but there is one 
case where a simple solution may be found. If the ini- 
tial conditions are thermal, as discussed above, then we 
should have a step function in c s in order to obtain scale- 
invariance (or require a 1). Let this step in c s (p) 
happen at p = p* and take c s = c_ (for p < p+) into 
c s = c + c_ (for p > p±). Integration of Eqn. (fT2"l) 
leads to w = c\ + (c 2 _ — c\)p/ p+ w c^_(l — p/p*) (for 
p > p*) showing that in this case we can ignore varia- 
tions in w while c s is changing, enforcing the conditions 
for our calculation while complying with (fT2"|) . Thus the 
calculations presented do not need to be modified if ther- 
mal initial conditions are used to obtain scale-invariance; 
in contrast with quantum initial conditions (and c s oc p), 
for which Eqn. (fT2"]l implies w w c^/3, so that a whole 
new calculation is warranted. 

We note that causality complaints [H, [3(| are bound 
to be model dependent. They don't affect bimetric and 
cuscuton models (see discussions in [13j, l3fj, |3l| ) , or more 
general K-essence models if they're seen as bimetric the- 
ories [3lJ. It may be that one generally needs to embed 
c s > 1 in bimetric VSL in order to prevent causality vio- 
lations; but the instantaneous change required in thermal 
scenarios could open up other possibilities. But it could 
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also be that VSL is not required to solve the causality 
paradoxes [32l |. 

4- Conclusions In summary we have revisited VSL sce- 
narios with reference to structure formation [i[ in the 
context of what we hope is a simpler framework: a vary- 
ing speed of sound. In VSL's initial formulation 0] 
the idea was to have increased symmetry as the Uni- 
verse cooled down, transitioning from a Galilean Uni- 
verse (with infinite speed of light) to the near Lorentzian 
Universe we see today. Obtaining a well defined formu- 
lation of such scenarios proved challenging, particularly 
under hard breaking of Lorentz invariance (for instance 
the issue of gauge choice became a physical one, and arbi- 
trariness ruled). Machian scenarios, where the constants 
of nature evolve alongwith the Universe, as a power of a 
were also considered 0, ^ is ironic that, as shown in 
this paper, thermal and quantum initial conditions pro- 
vide room for a phase transition [f| [|| and a Machian 
scenario [7j, respectively. 

How can we understand our results heuristically? In 
inflation a suitable heuristic can be obtained by not- 
ing that the curvature £ = v/z freezes outside the 
horizon and that inside the horizon v becomes a reg- 
ular Minkowski scalar field. The same is true in our 
scenario. In inflation one has to match free oscillat- 
ing modes of the form v ~ e lkv j\pZk inside the hori- 
zon with £ = F(k) cx fc~ 3 / 2 outside the horizon. The 
matching is done when kr\ ~ 1 and therefore requires 
z cx l/\rj\ for scale-invariance: a near deSitter back- 
ground. Here a similar argument can be made, but now 



inside the horizon modes have the well-known WKB form 
V ~ e lk J Cb I \/c s k, resulting from the variation in their 
frequency u = c s k. These have to be matched with 
C = F(k) cx k~ 3 / 2 outside the sound horizon, and the 
matching done when c s krj ~ 1. Thus we should have 
ayfc~ s cx resulting in the general condition In 
both cases study of the two extreme regimes (large and 
small) followed by suitable matching is enough to infer 
the final result. (Those conversant with Eqn. [7] — say, 
from inflation, or K-essence — will also quickly see why 
c s cx p leads to scale-invariance.) 

Are there any striking observational differences be- 
tween this mechanism and inflation? Deviations from 
n s = 1 can be easily obtained in this scenario, so we do 
not expect striking differences in terms of scalar fluctua- 
tions. However in the simple models wc have considered 
the horizon problem remains unsolved for gravitational 
waves and so we do not expect any tensor modes. This 
is to be contrasted with inflationary models with c s > 1 
where the tensor to scalar ratio is actually enhanced [ll| . 
This feature may be bypassed, say in more complex bi- 
metric theories, a matter we are currently investigating. 
Concrete implementations, such as the bimetric VSL the- 
ory, have also to be considered if one wants to address 
the other problems of Big Bang cosmology. 
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